Interband impact ionization and nonlinear absorption of terahertz radiations in 

semiconductor heterostructures* 
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We have theoretically investigated nonlinear free-carrier absorption of terahertz radiation in 
InAs/AlSb heterojunctions. By considering multiple photon process and conduction-valence in- 
terband impact ionization (II), we have determined the field and frequency dependent absorption 
rate. It is shown that (i) electron-disorder scatterings are important at low to intermediate field, 
and (ii) most importantly, the high field absorption is dominated by II processes. Our theory can 
satisfactorily explain a long standing experimental result on the nonlinear absorption in THz regime. 
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Many terahertz (THz) related physical phenomena 
HI HI [3 01 00 such as THz-radiation- induced dc cur- 
rent suppression multiphoton-assisted resonant tun- 
neling and multiphoton-assisted free-carrier absorp- 
tion [3j, have recently been reported. A few years ago, 
THz transmissions, reflections, and dc photoconductivity 
measurements were made to study multiphoton process 
and impact ionization (II) caused by THz radiations in 
InAs/AlSb heterojunctions (HJ). These experimental ad- 
vances motivated new theoretical studies 0, H El fiH El 
of THz physics. However, many interesting phenomena 
observed experimentally remained unexplained due to 
lack of well controlled theoretical models in the nonlinear 
regime. For example, the ac field dependent absorption 
coefficient increases with the field slowly at low field 
intensity, then increases rapidly at intermediate field in- 
tensities and levels off at high intensities. To date, this 
behavior is still not fully understood even within the pic- 
ture of multiphoton absorption mechanism. 

In this letter, we propose to explain this interest- 
ing absorption process for the first time by study- 
ing the multiphoton-assisted free-carrier absorption and 
electron- hole (e-h) generation in THz-driven InAs/AlSb 
HJ's. Most existing theoretical work emphasize the im- 
portance of the electron-impurity (e-i) or the electron- 
phonon (e-p) interaction 0, @ m nonlinear absorption. 
However, both these two scattering mechanisms lead to a 
decreasing absorption in the high field limit and thus un- 
able to explain the observed phenomena. The parameter 
that controls the absorption is the electro-optic coupling 
factor, r u = eE ac /(rau 2 ), here e is the carrier charge, 
to is the effective electron mass, E &c is the amplitude, 
and oj = 27r/ ac is the angular frequency with / ac the ra- 
diation frequency. In what follows we shall show that 
the dominant contribution to the nonlinear absorption 



in the intermediate to high field region is from the inter- 
band II processes. The calculated free-carrier absorption 
percentages of THz radiations in InAs/AlSb HJ are in 
excellent agreement with the experimental data 

Consider a semiconductor HJ with a two-dimensional 
(2D) energy- wavevector relation e s (kii), where s is the 
subband index. According to the balance-equation the- 
ory, when a uniform dc (or slowly varying) electric field 
Eo and a uniform sinusoidal radiation field, E(t) = 
Eo + E ac sin(27r/ ac i), are applied in the direction parallel 
to HJ interface, the transport state can be described by 
the following equations for the force, the energy, and the 
carrier number balance Ifl lift Il2r 
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where v is the average electron velocity, and JC is the 
inverse effective mass tensor. The electron sheet den- 
sity is N = 2 £ S(k|| / [(e,(k||) - n) /T e ] with T e the elec- 
tron temperature and fj, the electron chemical potential. 
f(x) = l/[exp(x) + l] is the Fermi distribution func- 
tion, g is the net electron-hole generation rate by bal- 
ancing II generation and Auger recombination. A e j, A ep , 
and An are the frictional accelerations respectively due 
to e-i scattering, e-p scattering, and interband II pro- 
cess. W ep and W\\ are the energy-loss rates respectively 
due to e-p scattering and interband II process. Si, S p 
and Su are the energy-gain rates of the electron system 
from the radiation field through the multiphoton process 
(n = ±I,±2,...) in association with e-i interaction, e-p 
interaction, and II process, respectively. The expression 
of Sji is as follow, 
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where T is the lattice temperature, q|| = (q x ,Qy) is the 2D 
phonon wave vector, U\ = k x v, L02 = wi + rnv 2 /2 + e£ + 
e S ,kn-q|| , w 3 = w 2 - ww, £ k|| = £ s , k|| - (1, ^ = - n h , 
£ k = Eg + kn/(2mh) is the hole dispersion with the 
effective hole mass, n(a^) = l/[exp(x) — 1] is the Bose 
function, and E g is the band gap. Il2(s', s, qn, fi) is the 
imaginary part of electron-electron correlation function 
0. J n (x) is the nth-order Bessel function. M], I s (q||) 
is the Fourier representation of the band-band Coulomb 
interaction matrix element |l3l| for the II and Auger pro- 
cesses in the 2D semiconductor system. The total energy- 
gain rate is defined by 



S = Si 



S p + Su. 



(5) 



It is useful to write S as a sum of all orders of n-photon 
contributions, S = Y^Li where S n is the total contri- 
bution from terms having index n and — n. We introduce 
the quantity a which is the ratio of the electromagnetic 
(EM) energy loss through the 2D sheet to the energy of 
the incident EM wave, a is a measure of the absorption 
of the radiation which can be written as, 



2NS 



(6) 



where c is the light speed in vacuum and k is the di- 
electric constant of the semiconductor. To see the role 
of individual multiphoton processes, we define a n = 
2NS n /{y/Ke cEl c ) and write a = Y^=i 

In the steady-state calculations for THz-driven 
InAs/AlSb HJ's, we consider the electron-acoustic- 
phonon scattering (via the deformation potential and the 
piezoelectric couplings) , electron-polar-optical-phonon 
scattering (via the Frohlich coupling), and elastic scat- 
tering both from the remote charged impurities and from 
the background impurities. For InAs, the Kane-type non- 
parabolic factor is 2.73 eV" 1 and E g = 0.22 eV. The 
three lowest electron subbands (e = 0, E\ — 35 meV, 
and e 2 = 200 meV) are included. Up to 10 hole sub- 
bands are taken into account. We set the sheet den- 
sity of InAs well No = 5 x 10 12 cm -2 , depletion layer 
charge density Ad cp = 5 x 10 10 cm -2 , background im- 
purity ni = 6.86 x 10 15 cm -3 , and remote impurities in 
AlSb barrier Nj — 1.53 x 10 11 cm -2 located at a distance 
of I = 10 nm from the HJ interface. In the whole paper, 
we set T = 300 K, and the dc field is assumed to be 
Eq = 1.5 V/m, same as that used in the experiments 0. 
The orders of the Bessel functions, n, in Eq. (J3J is set to 
be large enough so that the integral is convergent within 
the accuracy of 10 -4 . 



FIG. 1: (a) Electron-hole pair generation N/No and (b) elec- 
tron temperature T s /T vs i5 ac with and without II process in 
InAs/AlSb HJ at / ac = 0.64, 1, and 1.7 THz, respectively. 



In Fig. ^a) we show the e-h pair generation (solid 
lines) vs E ac with the II process at f ac — 0.64, 1, and 
1.7 THz, respectively. The dash line in Fig. da) shows 
no changes of the carrier number when the II process is 
excluded. The generation rate is essentially zero at low 
fields and high frequencies. If we define that the onset 
of the II process occurs at the threshold field E c (f ac ). 
The calculations indicate that E c (f ac )/f ac is nearly con- 
stant for all frequencies (see the inset of Fig. da), 
E c (f ac )/f ac ~ 3.59). The physical meaning of this lin- 
ear dependence is such that for II process to take place, 
times the photon energy fko must exceed a minimum 
value, r^huj = constant. At the high field iV depends 
on E ac exponentially. We shall show below that the ex- 
ponentially growing generation rate plays a crucial role 
in the total absorption percentage. The Fig. db) shows 
the effect of II processes on the electron temperature. 
While the mechanisms of the temperature increase with 
the coupling parameter r w is generally understood [lij . 
our results show a surprisingly large contribution to T e 
from the interband II processes. The interband II pro- 
cesses account for about 50 percent of T e . It is interesting 
to note that II processes lead to a net increase of T e . In 
many cases, II processes act as a cooling mechanism. For 
example, for a 3D system under a dc field, T e with II pro- 
cess is lower than that without II process. In this case 
the electron heating is only due to electron accelerating 
and scattering under the dc field. For the present system 
under an ac field, the II processes still act as a cooling 
mechanism. However, large number of electrons created 
in II processes immediately causes a rapid increase in 
photon absorption due to e-i, e-p, and more importantly, 
e-h scattering. This absorption is further enhanced by 
nonlinear multiphoton absorption under an ac field. This 
subsequent electron scattering and absorption of photon 
energy lead to an increase of T e which exceeds the initial 
reduction of T e due to II processes. As a result, T e is 
higher with II processes. 

In Fig. |2fa) we show the calculated total absorption 



3 



percentage a vs E ac at / ac = 0.64, 1, and 1.7 THz, re- 
spectively. The solid circles are the experimental results 
for / ac = 0.64 THz from Ref. 3, where a = 1 - T r - R 
with T r the transmission and R the reflection. Excel- 
lent agreement is obtained between the calculated results 
and the experimental data. It's seen from Fig. [2Ja) 
that, lower frequency and/or higher radiation intensity 
lead to more absorption. For / ac = 0.64 THz, the ab- 
sorption rate increases slowly from low electric field to 
about E ac — 3 kV/cm, then goes fast before reaching its 
saturation value. The significance of II process can be 
clearly seen in the measurements and in the current the- 
ory. Without the II process, the multiphoton absorption 
is due to e-i scattering and e-p scattering. The absorption 
rate due to these mechanism only increases very slowly at 
low E ac . In Fig. Efb) we show the phonon-induced (S* p ), 
Il-induced (Su), and the total energy gains (S p + Su) 
of the electron system from the radiation fields. At in- 
termediate to high E ac , both e-i and e-p scatterings are 
actually suppressed by the strong electro-optic coupling. 
As a result one observes a rapidly decreasing absorption 
in the high field limit. As we have seen in Fig. 2] the 
effect of strong field on the II process is quite opposite 
to that of electron disorder scattering because the e-h 
pair generation is a fundamentally different process than 
the scattering process. Scattering depends on the abil- 
ity of electrons to simultaneously exchange momentum 
and energy with the scatters and photon. Strong cou- 
pling with photons reduce the probability of momentum 
exchange with the impurities or phonons. This momen- 
tum exchange is not required in an e-h pair generation 
process. Annihilation of multiple photon in a genera- 
tion process leads to a strong total absorption percent- 
age. The very rapid increase of II process at high field 
[see Fig. ^a)] is largely offset by the fast decreasing 
of impurity and phonon scattering. The resulting to- 
tal absorption percentage increase modestly in the high 
field limit. The curve at / ac =0.64 THz can be directly 
compared to experimental result. At this frequency, II 
mechanism begins to make contribution to the total ab- 
sorption after about E ac = 3 kV/cm. The carrier number 
increases with increasing E ac . For / ac = 0.64 THz, when 
E ac = 6.5 kV/cm, the total electron number N increases 
up to about six times the initial number Nq due to the 
II process. The agreement between the present theory 
and the experiment is excellent. In the inset of Fig. Gib) 
we show the n-photon absorption percentage a n vs the 
order n of the Bessel function at / ac = 1.7 THz. E ac is 
changed from 2 to 6 kV/cm with a step of 0.5 kV/cm. 
For each E ac there is a maximal of a n , and the posi- 
tion of the peaks shifts to higher order n with increasing 
E ac . It means that higher order multiphoton processes 
play more important role on the total absorption with 
increasing E ac . Generally, it can be seen that the role 
of THz radiations with larger E ac and lower / ac on car- 
rier transport increases with increasing the electric-optic 



coupling factor r u . 

In Fig.|3 we show the total absorption percentage a 
vs the wavelength A at E ac — 2 kV/cm (triangles), 5.8 
kV/cm (circles), and 9 kV/cm (squares), respectively. 
The absorption percentage a increases with increasing 
E ac and/or increasing A. The relation between a and 
A can be fitted by the function (see the lines in Fig. 
EJ: a(X) = a X b . When E ac = 2 kV/cm the fitted 
function is a(A) = 1.70 x 10~ 6 A 2 , which is the classi- 
cal absorption-wavelength relation. When E ac increases, 
however, the present nonlinear theory predicts that the 
square relation will be broken. For E ac — 5.8 kV/cm we 
have (a , b) = (7.19 x 10~ 8 , 2.6), and for E ac = 9 kV/cm 
we have (an, b) = (3.31 x 10~ 8 , 2.8). The physical reason 
is the exponential growth of the e-h pair and the more im- 
portance of the multiphoton process with increasing r u . 
In the inset of Fig. |3| we show the first five order multi- 
photon contribution a„ (n = 1,2, ...,5) vs A at E ac = 2 
kV/cm. It can be seen that, for each order of multipho- 
ton channels there is a critical wavelength A p at which 
the multiphoton process makes maximal contribution to 
the total absorption. With increasing n, the position of 
A p shifts to longer wavelength. The physical origin of the 



FIG. 2: (a) Absorption percentage a vs E ac at / ac =0.64, 1, 
and 1.7 THz, respectively. The circles are the experimental 
results from Ref. 3. (b) The corresponding phonon-induced 
(S p ), Il-induced (Su), and total energy gains (S p + Su) vs 
-Eac- The inset shows n-photon absorption percentage a n vs 
the order n of the Bessel function at / ac = 1.7 THz. 
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FIG. 3: Calculated total absorption percentage a vs wave- 
length A at _E ac = 2 kV/cm (triangles), 5.8 kV/cm (circles), 
and 9 kV/cm (squares), respectively. They are fitted by the 
function (lines): ot(\) — aoA 6 . The inset shows the first five 
order multiphoton contribution at E ac — 2 kV/cm. 



critical wavelength at a given photon number is that, for 
a given n, the strength of the scattering matrix of the 
n-photon assisted process decreases almost like u>~ 4n for 
large u>, while the energy absorpted from the radiation 
field during each n-photon process is proportional to nuj 
[see Eq. Q]. There must be a critical wavelength X p at 
which the energy absorption reaches a maximum. 

Several physical mechanisms employed in our theoret- 
ical model and numerical calculation need to be further 
discussed. (1) We have neglected the absorption of bulk 
electrons in the substrate because the total absorption 
is dominated by the 2D channel. The Coulomb inter- 
action matrix for the 2D channel is much stronger than 
the corresponding interaction matrix of 3D system ex- 
cept in the small wavevector regime. On the other hand 
the electron-laser coupling qy • r w [sec Eq. (@J] is only im- 
portant at large qu regime. As a result, the THz-induced 
II in the bulk plays a negligible role in the present sys- 
tem. (2) Although a large number of holes are being 
created in II processes, the hole contribution to the total 
absorption is ignored in our model. The hole effective 
mass (mh = 0.94mo) in our system is about 25 times 
heavier than that of the electron (m = 0.038mo), here 
mo is the mass of the free electron. From Eq. (4), one 
can see that Su is mainly determined by J^(qii • For 
a small r u , i.e., larger mass m or higher frequency oj, 
we have J 2 (q|| • r w ) oc m~ 2 l™ for \n\ > 0. For n = 0, 
we can directly see from Eq. (4) and the inset of Fig. 



|5Jb) that the term does not contribute to the total ab- 
sorption. When the electron mass is replaced by the 
hole mass in J 2 (q|| ■ r w ) (for |n| > 0), this quantity 
will be reduced by a factor of around 25 -2 = 0.0016. 
This makes the hole contribution at least two to three 
orders of magnitude smaller than that of the electrons. 
(3) We have used a strict 2D model for charged carriers 
and neglected the effect of real space transfer (RST) to 
bulk electrons. For InAs/AlSb HJ svstem, the confin- 
ing potential is AE C — 1350 meV |l5j. This is about 
15666 K. The maximum T e in our calculation is around 
20 times the lattice temperature (about 6000 K). Since 
the case / ac = 0.64 THz represents our central result, 
we estimate the RST effect. At E ac = 6 kV/cm, T e 
is about 12T or 3600 K, which is less than a quarter 
of the confining potential. If the RST is due to over- 
the-barrier activation, the probability of this transfer is 
~ exp(-AE c /T e ) « 0.0129, which is less than 1.5 per- 
cent. On the other hand, in terms of photon energy of 
1 THz (about 4 meV), the barrier height is about 330 
times of photon energy. This makes the RST by gaining 
enough energy from the THz field also negligible. Ac- 
tually, it is clear from the inset of Fig. Efb) that the 
probability of absorbing more than 30 photons is practi- 
cally zero. Therefore we neglect the RST effect. 

In conclusion, we have shown that II process plays an 
important role in the free-carrier absorption of THz radi- 
ations in InAs/AlSb HJs. Our results can explain a long 
standing experimental result on the free carrier absorp- 
tion in the THz regime. 
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